Search for neutrino oscillations at a research reactor 
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A project of search for oscillations of reactor neutrinos on a short, 5-15 meters, baseline with liquid 
scintillation position-sensitive detector POSEIDON is described. The oscillations of the electron 
antineutrinos into sterile neutrinos will be searched at 100 MW research reactor with the size of the 
active zone of about 0.5 m. The values of the oscillation parameters available for research are in the 
region Sm 2 = (0.3 - 6) eV 2 and Sin 2 (2<9) > 0.01. 
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I. INTRODUCTION 



At present there is compelling evidence for the exis- 
tence of sterile neutrino with the mass splitting Sm 2 ~ 
1 eV 2 relative to three known mass states and mixing pa- 
rameter Sin 2 (29) ^0.1. Since the decay width of the 
Z-boson indicates that the number of light neutrino flavor 
states does not exceed three, a new sterile neutrino has 
no weak interaction. The discovery of the new fourth (or 
even fifth) neutrinos would be beyond the scope of SM 
and would mean the existence of new physics unknown to 
us. The possible presence of one or two sterile neutrino 
states points towards non-standard neutrino physics be- 
yond the Standard Model (SM). 

The schemes with a fourth massive neutrino with a 
mass of ~ 1 eV 2 have attracted attention since the LSND 
group announced an excess of electron antineutrinos in 
muon antineutrino beam [l| . This evidence has been con- 
firmed by the MiniBoone experiment in the antineutrino 
and the neutrino sector Q. The data from the KAR- 
MEN detector exclude the region of large values dm 2 > 
(~ 5 eV 2 ) i. 

The other indications of short oscillations with 
8m 2 ~ 1 eV 2 were obtained by radiochemical Ga-Ge- 
detectors SAGE and GALLEX. The measurements with 
the artificial neutrino sources 51 Cr and 37 Ar 0, Q in- 
dicate a deficit of electron neutrinos and it may imply 
neutrino oscillations. 

Recently, the new calculations of the spectrum of re- 
actor neutrinos lead to the value of the neutrino flux by 
about 3% greater than which is consistent with the ex- 
periments H, 0|- As a result, the data of the previous 
reactor experiments can be interpreted as the observa- 
tion of a flux deficit, and this, in turn, can be explained 
by additional sterile neutrinos with the masses at the eV 
scale 

At last, the current data from cosmology and Big 
Bang nucleosynthesis favor an additional relativistic de- 
gree of freedom which can be interpreted as sterile neu- 
trino species 

Some experiments to search for oscillations with 
Sm 2 ~ 1 eV 2 were proposed at reactors and at accel- 
erators and with artificial neutrino sources. The value 
8m 2 ~ 1 eV 2 corresponds to the ratio of oscilla- 



tion length to neutrino energy L/E^ w 1 m/MeV = 1 
km/GeV. 

The Borexino collaboration plan to perform measure- 
ments with artificial sources of neutrino ( 51 Cr, 37 Ar) and 
antineutrino ( 144 Ce, 90 Y+ 90 Sr) which makes it possible 
to check further CPT invariance f^-fllj]. The Borexino 
detector registers the scattering of neutrinos by electrons 
and the reaction of inverse beta-decay for antineutrinos 

The possibility of observing neutrino oscillations with 
51 Cr source, through the allocation of space zones of Ge 
production, is considered for the detector SAGE [l4j . 

The experiment with the sources of 144 Ce and 106 Ru, 
which can be placed in the center of a large volume scintil- 
lation detector such as SNO, BOREXINO or KamLAND 
is proposed in [15j. Detectors Nucifer (Saclay/IN2P3) 
[L6J and DANSS (JINR / ITEP) ^ which are created to 
control the nuclear fuel, plan to use to search for oscilla- 
tions at short distances. The experiments at accelerators 
with L/E„ ~ 1 km/GeV are also considered fl8j . 



II. 



DETECTION OF ANTINEUTRINO 



The reaction of inverse beta decay v + p — > e + + n is 
used for the registration of the antineutrino. Two con- 
secutive events from the positron and neutron capture in 
a short time window, which is determined by the lifetime 
of a neutron in the scintillator (^25 ms with 0.1% Gd 
mass), can reliably identify this reaction. The thresh- 
old of the reaction is 1.8 MeV and the positron energy is 
uniquely related to the antineutrino energy. At a distance 
of 5 meters from 100 MW reactor the expected number 
of (vp, ne + ) events is 3 / (kg day) for a target containing 
6xl0 25 protons / kg (1 kg C3„H4„ scintillator). 

The length of the oscillations for the two neutrinos 
mixing is: L [m] = 2.5 E„ [MeV] / Sm 2 [eV 2 ], and the 
depth of the oscillations is determined by the angle of 
mixing Sin 2 (28) . The probability of observing the elec- 
tron antineutrino (P oc ) passing the path R is described 
by simple expression P oe = 1 - Sin 2 (29) Sin 2 (R/L). For 
a non-point source, as well as for a non mono-energetic 
neutrinos, P ee does not depend on R for R 3> L: P oc = 1 - 
0.5 Sin 2 (29). In this case, the oscillations appear through 
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FIG. 1. The value of o-(vp, ne )xP ee depending on the dis- 
tance R for the neutrino detection threshold 1.8 MeV. 



FIG. 2. The value of cr(i/p, ne + ) xP ee depending on R for the 
different ranges of the neutrino energy. 



a decrease in the detected neutrino flux as compared with 
the flux calculated for the absence of oscillations. 

We propose to look for the oscillation waves in the 
counting rate of antineutrinos of different energies asso- 
ciated with the transitions of the electron antineutrino 
in a sterile one v e — > v s and v s — > v e , directly 
on the size of the detector. Two conditions are required 
to perform this. The size of the neutrino source must be 
smaller than the length of the oscillations. Since the re- 
actor neutrinos have a continuous spectrum, the energy 
resolution of the detector should be sufficient to observe 
changes in the shape of the neutrino spectrum. The os- 
cillation curves on the length of the detector are shown 
in Figures 1 and 2 for the values 8m 2 = 1.5 eV 2 and the 
angle of mixing Sin 2 (2#) = 0.15. These values are close 
to the ones obtained from the joint fit of the reactor, gal- 
lium an MiniBooNE data (8m 2 > 1.5 eV 2 , Sin 2 (26>) = 
0.14 ± 0.08) 0. 

In the case of registration all neutrinos, with the en- 
ergies above (vp, ne + ) reaction threshold, the deviation 
from the law 1/R 2 do not exceed 1.8%, if the neutrinos 
are detected at a distance of more than 6 m (Fig. 1). But 
for the antineutrino with the energy in the range (3 - 4) 
MeV, the change is 9% over the length of the detector 
1.5 m (Fig. 2). To observe these oscillations waves, the 
energy resolution of the detector should be better than 
0.5 MeV. 



The most intense sources of neutrinos are nuclear 
power plants with thermal power of ~ 3 GW, which pro- 
vides the neutrino flux 10 13 cm _2 s _1 at the distance of 
20 meters. The characteristic size of the zone is about 5 
meters that makes it impossible to observe the oscillation 
waves if 8m 2 > 1 eV 2 . The zone of research reactors is 
much smaller and is about 0.5 meters. Thus, the zone of 
the high-flux reactor PIK in Gatchina has the height of 
50 cm and the diameter of 39 cm [l^| . The zone of the re- 
actor SM-3 in Dimitrovgrad is even smaller - 42 x 42 x 35 
cm [H. 



III. 



THE POSITION-SENSITIVE DETECTOR 



The proposed detector with the dimensions of central 
part 1.5 x 1.0 x 1.0 m 3 is the liquid scintillator (LS) detec- 
tor. The scintillation signals are detected by 72 photo- 
multipliers PMT125 with the flat photocathode 150 mm 
in diameter. PMTs are placed on two opposite yz-planes 
(if the reactor active zone lies on x-axis) of the detec- 
tor and they are separated from the central volume by 
30 cm buffer. The buffer (no scintillation) provides the 
spectrometric properties of the detector and protects the 
detector from radioactivity of PMTs. The central de- 
tector is surrounded by 4 additional 1.5 x 1.2 x 0.15 m 3 
detectors, each of which is viewed from the both sides 
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by 14 PMTs. The sections perform the role of the active 
shielding against cosmic rays and fast neutrons. They 
also detect 511 keV gammas which escape from the cen- 
tral detector after the positron annihilation. The scintil- 
lators based on the linear alkylbenzene [Il|,[22[ (LAB, T 
= (130-150) C°) or PXE [H (Ci 2 Hi 8 , T = 167 C°) are 
preferable because of a high flash point. The properties 
of PXE were studied in details on a prototype Borex- 
ino detector (23|. LAB is used by Daya Bay 24 1 and 
Reno [25| reactor experiments. Both solvents allow the 
production of the Gd-containing scintillators. 

The energy and coordinates of an event are recon- 
structed from the charge registered by PMTs. The Monte 
Carlo simulations based on GEANT4 were performed to 
find the main parameters of the detector. Fig. 3 shows 
the detector response functions for 1 MeV (kinetic en- 
ergy) positrons. The positrons were uniformly simulated 
inside the central detector. A light yield equal to 10 4 
photons/MeV is a real value for the scintillators based 
on the LAB or PXE. For the quantum efficiency of the 
photocathode PMT125 equal to 77 = 0.1, about 200 the 
photoelectrons (p.e.) at 1 MeV energy release are regis- 
tered, which provides a statistical component of the en- 
ergy resolution of a = 7% / E 1 / 2 where the energy E is 
expressed in MeV units. The 511 keV escape peak is sig- 
nificantly suppressed for the case of additional sections 
of a thickness of 15 cm (Fig. 3, curve 2). 

To restore the coordinates of the event, the values of 
the charges (Qi, i = 1 .. 72), registered with 72 PMTs are 
used. The reconstruction program is based on a simple 
ratio of the charges collecting by PMTs from front and 
back sides of the detector. It provides a spatial resolution 
of 15 cm (la) for the kinetic energy of positrons 1 MeV 
(Fig. 4). The reconstruction of the x-coordinate (in the 
direction of the reactor core) with an accuracy of bet- 
ter than 15 cm allows observing oscillation waves on the 
length of the detector. At the same time this allows to 
improve the energy resolution due to the space position 
and to compensate for the loss of the detection efficiency 
in the edges of the detector. 

The spatial resolution for the neutron (a = 25 cm) 
is significantly worse than for the positron, however, it 
can be used as an additional criterion for the selection 
of inverse /3-decay events. In addition, time signals from 
PMTs allow us to distinguish the signals from electrons, 
cv-particles and recoil protons due to different photons 
time responses. 

The background level is the main problem of the ex- 
periment carried out at the Earth's surface and at a min- 
imum distance from the reactor. The lead shield with the 
thickness of 200 g/cm2 suppresses the external gamma- 
ray activity almost completely. The copper or iron shield 
is undesirable because of the high-energy gamma rays 
from the capture of thermal neutrons. The hydrogen- 
consisting shield such as borated polyethylene reduces 
the flux of fast and thermal neutrons inside the detector. 
All construction materials are tested for radioactive pu- 
rity in the low-background setup with HPGe-detector. 




FIG. 3. The detector response functions for 1 MeV positrons: 
1 - in the central volume ; 2 - in the whole volume, 3 - in 15 
cm outer layer. 



The cosmic rays background is reduced by the active 
shield constructed from plastic scintillators. The rela- 
tively small size of the main detector allows to use 2 
layers of the active shield before and after the lead shield 
which is the generator of neutrons. 

PMTs operate in a mode close to one-electron regime. 
The amplifier has a short restore time after pulses corre- 
sponding to a muon passing through the LS. The count 
rate for a single PMT set about 100 Hz. The trigger ap- 
pears if more than N PMTs fire in the time window of ~ 
20 ns. E.g. N = 20 corresponds to the energy threshold 
of 100 keV. The states of all PMTs are recorded in the 
time range 50 ms before and after the trigger. The charge 
and time signals are used to reconstruct the position and 
energy of the event. For the trigger rate equal to 10 s _1 , 
about 1 GB of data are stored per day. The neutrino 
events are selected in the off-line mode. 



IV. THE MAIN SOURCES OF THE 
BACKGROUND 

The main difficulty of this experiment is that it is car- 
ried out almost on the surface of the Earth. The sup- 
pression of the muons and the nuclear-active component 
of cosmic rays is with only a few (less than 10) meters of 
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FIG. 4. The restore of x-coordinate for the events simulated 
in the center of the detector. 



water equivalent. The detector background is the main 
problem of the experiment. A relatively small neutrino 
flux will be detected during reasonable time only if the 
background is comparable to the effect. 

The main source of the background not correlated with 
the reactor is induced by random coincidences within a 
time window r when signals from positron and neutron 
are selected. The number of random coincidences is de- 
termined by the value N1N2T where Ni and N2 are the 
count rate of in the ranges where positrons and neutrons 
are registered. The value r = 25 ms obtained for a 0.1% 
Gd concentration can not be reduced significantly. As a 
result, the value N1N2 must to be less than 300 s~ 2 , in 
order to the background of random coincidences does not 
exceed the number of inverse beta-decays. To decrease 
N1N2 it is needed to increase the registration thresholds 
for positron and neutron or to sectionalize the detector. 

The background due to cosmic rays is suppressed by 
an active shield. Since the count rate of veto is at least 
10 3 s" 1 , the length of veto can not exceed 100 ms to 
avoid a large dead time. As a result, the neutrons pro- 
duced by cosmic rays can give a signal behind the veto 
pulse. Another source of the background induced by 
long-lived light nuclei ( 8 Hc, 9 Li, etc) which decay with 
neutron emission. The rate of the production of such 
nuclei by muons in the Earth's surface is ~ 150 events 
/ (t day). The production of long-lived light nuclei by 
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FIG. 5. The count rate multiplied by the value (R [m] / 6 m) 
2 . The errors correspond to 100 days of data taking. Dotted 
line is no-oscillation case. 

the nuclear-active component will be measured with the 
detector prototype. 

The natural radioactivity of the detector materials and 
the first layer of the passive shield determines the regis- 
tration threshold for positrons. The 40 K in a glass PMT 
is the main source of the background. Due to 40 K the 
threshold for positron will be set above 1.5 MeV that cor- 
responds to the minimum neutrino energy E u = 2.3 MeV. 
Due to the radioactivity of the lead shield and other inner 
components, the threshold may be increased beyond the 
natural gamma activity of 2.6 MeV (E„ = 3.4 MeV). The 
total energy of the gammas produced by neutron capture 
by Gd nuclei is 8 MeV. The MC simulations show that 
the threshold of 2.6 MeV for the delayed signal corre- 
sponds to the detection efficiency of 80%. 

To suppress the background of fast neutrons, corre- 
lated with the reactor power the difference in time re- 
sponse of the scintillator for electrons and the recoil pro- 
tons is used (/3/p discrimination). This discrimination 
suppresses the background of successive Bi-Po (3 - and 
a-decay in uranium and thorium families. 

V. SENSITIVITY OF EXPERIMENT 

We defined the sensitivity of the experiment for the 
neutrino energy range (3 - 4) MeV where the count rate 
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is 900 events per day if the front wall of the detector 
is located at a distance of 6 m from the reactor core. 
Fig. 5 shows the count rate of the reaction of inverse 
beta decay depending on the distance from reactor core. 
The calculations take into account the finite size of the 
active zone and the energy and spatial resolution of the 
detector. If Sin 2 (2(9) = 0.15 the difference in the number 
of counts for front (R = (6-6.75) m) and back (R = (6.75- 
7.5) m) volume of the detector is 4.5%. By assuming 
that the background is equal to the effect and detection 
efficiency close to 100% about 200 days of taking data 
are needed to measure Sin 2 (2$) with an accuracy of 0.04. 
The number of events for a neutrino detected in the range 
3-4 MeV is only 28 % of the full neutrino events. 

The proposed scheme of the experiment is realized only 
if the oscillation length lies in the range (1.5 - 30) m. The 
values of the parameters of neutrino oscillations, which 



can be studied using the deviations from 1/R 2 , lie in the 
region Sm 2 = (0.3 - 6) eV 2 and Sin 2 (26>) > 0.01. The 
sensitivity and reliability of the experiment increase by 
using two identical detectors Dl and D2 at two distances 
of 6 and 9 meters (Fig. 6). 

VI. CONCLUSION 

The project with a liquid-scintillation position- 
sensitive detector (POSEIDON -position-sensitive detec- 
tor of neutrino) to search for oscillations of reactor neutri- 
nos at short (5 - 15) meter base is proposed. The project 
aimed to search for oscillations of the electron antineu- 
trinos into sterile ones at 100 MW research reactor with 
0.5 m active zone. The oscillation parameters which can 
be investigated are in the range 5m 2 = (0.3 - 6) eV 2 and 
Sin 2 (26>) > 0.01. 
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